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ABSTRACT 

Understanding the properties of Pop III stars is prerequisite to elucidating the nature of primeval 
galaxies, the chemical enrichment and reionization of the early IGM, and the origin of supermassive 
black holes. While the primordial IMF remains unknown, recent evidence from numerical simulations 
and stellar archaeology suggests that some Pop III stars may have had lower masses than previously 
thought, 15 - 50 Mq in addition to 50 - 500 Mq. The detection of Pop III supernovae by JWST, 
WFIRST or the TMT could directly probe the primordial IMF for the first time. We present nu- 
merical simulations of 15 - 40 Mq Pop III core-collapse SNe done with the Los Alamos radiation 
hydrodynamics code RAGE. We find that they will be visible in the earliest galaxies out to z ~ 10 - 
15, tracing their star formation rates and in some cases revealing their positions on the sky. Since the 
central engines of Pop III and solar-metallicity core-collapse SNe are quite similar, future detection of 
any Type II supernovae by next-generation NIR instruments will in general be limited to this epoch. 

Subject headings: early universe - galaxies: high-redshift - stars: early-type - supernovae: general - 
radiative transfer - hydrodynamics - shocks 



1. INTRODUCTION 

The first stars are crucial to the formation of primeval 
galaxies, the chemical enrichment of the early IGM, 
the initial stages of cosmological reionization, and the 
origin of supermassive black holes. Unfortunately, be- 
cause they lie beyond the reach of current ground and 
space based instruments there are not yet any obser- 
vational constraints on their properties. The early nu- 
merical simulations of Pop III stars suggest that they 
form in isolation, one per halo, in 10 5 - 10 6 Mq dark 
matter halo s at z ~ 20 - 30 and are very massive, 100 

500 M m (IBromm et all H99l lAbel etaTJl20M [20021: 



iBromm et alj|2002t iNakamura fc Umemura 120011 ) . Usu- 
ally, such stars drive most of the baryons from their 
halos in strong ionized flows that create diffuse H II 
regions with n ~ 0. 1 - 1 cm" 3 (IWhalen et all [2001 
Kitayama et all 120041: lAlvarez et al.l 120061: lAbel et al l 



2007tlWise fe AbeW2008ft . Pop III stars are not thought 



to lose much mass over their lifetimes because there 
are no line- driven winds in their m e tal-free atmospheres 
(IKudritzkil [20001: IVink etHl l200lt IBaraffe et alTl200l"l: 
lEkstrom etaiTl2008t) . 

There is growing evidence that some Pop III stars may 
be less massive than previously thought. Recent, more 
extensive ensembles of numerical simulations have found 
many halos with central collapse rates consistent with 20 
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- 60 M q for the final mass of the star (jO'Shea fe Normanl 
12007ft and th at a fraction of th e halos form binaries in this 
mass range ([Turk et all 12009ft . New simulations of Pop 
III protostellar accretion disks suggest that they were 
prone to fragmentation i nto as many as a dozen smaller 
stars (IStacv et all 120101: IClark et all [20Tlt ISmith et all 
120111: iGreif et all 120111 12012ft . Preliminary models also 
suggest that ionizing UV breakout from primordial pro- 
tostellar disks limits the masses of s ome Pop III stars 
to ~ 40 M Q (IHosokawa et all 120111: IStacv et all 12012ft 
(alth ough also see iTan fe McKed 120041 : IMcKee fe Tanl 
[2008ft . 

While these new models are important steps forward 
they cannot predict the final masses of primordial stars. 
For example, the fragmentation of Pop III protostellar 
disks has only been followed out to a few centuries, far 
short of the time required to build up a massive star, 
and many of the fragments in these models are later 
found to merge with the central object instead of be- 
coming stars themselves. Thus, it is not clear if a small 
swarm of less massive stars forms or a single, very mas- 
sive star is created, albeit through protracted clumpy 
accretion. Clumpy accretion may also keep the proto- 
star puffed up and cool, allowing it to reach much higher 
masses before evapora ting its accretion disk than in the 
IHosokawa et al.l (|2011ft models. Because no simulation 
realistically bridges the gap between the formation and 
fragmentation of the accretion disk and its destruction 
up to a Myr later, they cannot yet constrain the Pop III 
initial mass function (IMF). 

There have been attempts to determine the masses of 
the first stars by comparing the cumulative nuclcosyn- 
thetic yield of their supernovae (SNe) to the fossil abun- 
dance record, the chemical abundances measured in an- 
cient, dim extrem ely metal-poor (EMP) stars out in the 
Gala ctic halo (e.g. IBeers fe Christliebll2005l: iFrebel et ail 
12005ft . Pop III stars from 15-40 Mq die in core col- 
lapse (CC) SN explosions and 140 - 260 Mq stars ex- 
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plode as far more ener getic pair instability (PI) SNe 
(|Heger fc Wooslevll2002D (although the lower mass limit 
for PI SNe has recent ly dropped down to 65 Mp, fo r 
rotating progenitors - iChatzopoulos fc Wheeler! I2012h . 
Uoggerst et all (|2010t ) (hereafter JET10) found that the 
chemical yields of 15 - 40 Mq Pop III SNe are a good 
match to th e abundances observed in a sampl e of ~ 130 
EMP stars (jCavrel et al.ll200l lLai et alll2008T) . and that 
such stars may have been responsible for much of the 
chemical enrichment of the primeval IGM. This, together 
with the failure to detect the distinctive odd-even nucle- 
osynthetic imprint of PI SNe in EMP stars to date, has 
led some to conclude that Pop III stars could not have 
been very massive. However, second-generation stars 
may have been enriched by the ashes of the first PI SNe 
to metallicities well abov e those targeted by surveys thus 
far (jKarlsson et alJl2008f > , so such conclusions are prema- 
ture. In sum, while stellar archaeology has revealed some 
insights into the first stellar populations it has not yet 
placed firm constraints on their masses. 

Observations of primordial SNe could directly probe 
the Pop III IMF for the first time because in princi- 
ple they can be detected at great distances and distin- 
guish between low mass and high mass progenitors. Until 
now, calculations of Pop III SN light curves and spec- 
tra have been confined to PI SNe to determine if they 
can be detected by futu re SN la missions (4 < z < 6, 
IScannapieco et all 120051) or by the Ja mes Webb Space 
Telescope (JWST. [Gardner et alll2006l) during the era of 
reion ization (6 < z < 15, Kasen et all 1201 ltlPan et al 



2012]) or at z ~ 30 (IWise fc AbeU 12005b iHummeT et a 



20121 ). Simulations of Pop III PI SN light curves have also 



been applied t o PI candidates in th e local universe suc h 
as SN 2007bi (|Gal-Yam et all 120(1 : lYoung et alJl2010h . 
Most recently, radiation hydrodynamical models of Pop 
HI PI SN light curves and spectra by iWhalen et al.l 
(|2012d lbT) showed for the first time that these explosions 
will be visible out to z ~ 30 with JWST and z ~ 15 
- 20 in all-sky near infrared (NIR) surveys such as the 
Wide-Field Infrared Survey Telescope ( WFIRST) or the 
Wide-Field Imaging Surveyor for High-Redshift ( WISH). 
However, no calculation to date has addressed detection 
thresholds in redshift for Pop III CC SNe, whose numbers 
may be greater than those of PI SNe at early epochs. 

We present radiation hydrodynamical simulations of 
light curves and spectra for 15-40 Mq Pop III explo- 
sions done with the Los Alamos RAGE and SPECTRUM 
codes. In § 2 we describe our grid of explosion models and 
how we post process them obtain light curves and spec- 
tra. In § 3 we examine blast profiles and spectra for Pop 
III CC SNe in detail, and in § 4 we compute their near 
infrared (NIR) light curves and detection thresholds in 
redshift. In § 5 we conclude. 

2. NUMERICAL ALGORITHM 

We adopt the models in JET 10 for our grid of light 
curve simulations because they span the range of progen- 
itor masses, structures and explosion energies expected 
for such stars and because their elemental yields are a 
good match to those in the fossil abundance record. We 
calculate light curves and spectra in four stages. First, 
15-40 Mq Pop III stars are evolved through all stages of 
stable nuclear burning and then exploded in the Kepler 
code. Each explosion is then mapped onto a 2D adaptive 



mesh refinement (AMR) grid in the CASTRO code and 
evolved until just before the shock reaches the surface of 
the star to capture internal mixing prior to breakout. We 
then spherically average our CASTRO profiles onto ID 
AMR grids in RAGE and evolve them out to 4 months, 
after which the ejecta dims below observability. Finally, 
we post process our RAGE profiles with the SPECTRUM 
code to construct light curves and spectra with detailed 
atomic opacity sets. 

2.1. Kepler 

We use the ID Lagrangian stellar evo lution code Ke- 
pler (jWeaver et al.lll978t IWooslev et all 120021) to evolve 
the progenitors from the beginning of the main sequence 
to the onset of collapse of their iron cores. The explosions 
are then artificially triggered by a piston at a constant 
Lagrangian mass coordinate that advances though the 
star with a specified radial history. The SNe are evolved 
through all nuclear burning (~100 s after the explosion) 
and then the simulations are halted, well before the shock 
exits the He shell and any reverse shocks form in which 
instabilities might develop. We follow energy production 
with a 19-isotope network until oxygen is depleted from 
the core and with a 128-isotope quasi-equilibrium net- 
work thereafter. Stellar rotation is approximated by a 
semi-convective mixing parameter. 

As in JET10, we consider two progenitor metallicities 
(Z = and 10 -4 Zq denoted by z and u, respectively), 
three explosion energies (E ex = 0.6, 1.2, and 2.4 B des- 
ignated by B, D, and G, respectively, where 1 B = 1 
Bethe = 10 51 erg), and three progenitor masses (15, 25, 
and 40 Mq), a total of 18 models. Thus, model zl5G 
is the 2.4 B explosion of a zero-metallicity 15 Mq Pop 
III star. We use only one rotation rate from that study 
(R = 5, where R is defined in § 3 of JET10) because it 
was found that the rate of rotation had little effect on 
the degree of mixing during the explosion. In contrast to 
the u- and z-series 150 - 250 Mq stars in IWhalen et al.l 
(20121)), u-series 15 - 40 Mq Pop III stars die as compact 
blue giants and z-series stars die as red supcrgiants. 

We refer our readers to § 2 of JET 10 for further details 
on our Kepler models, such as our numerical mesh, ex- 
plosion triggers, and use of semi-convective mixing as a 
proxy for rotation. 

2.2. CASTRO 

As in JET10, we map our Kepler explosion profiles 
onto a 2D axisymmet ric grid in the new CASTRO code 
(jAlmgren et alJ 120101 ) and evolve the shock up to the 
edge of the star in mass coordinate. CASTRO is a mul- 
tidimensional Eulcrian AMR code with a high order, un- 
split Godunov hydrodynamics scheme. We follow the 
same 27 elements from hydrogen to zinc as in JET10 
and include the gravity of the compact remnant, which 
JET10 and others have shown to be crucial to both fall- 
back and mixing in Pop III CC explosions (jZhang et al.1 
120081 ). Radiation transport is not necessary in this stage 
of our calculations because photon mean free paths in the 
star prior to breakout are so short that they are simply 
advected through the star by the fluid flow. We do ac- 
count for the contribution of photons to pressure in the 
equation of state (EOS). Our models also include energy 
deposition due to radioactive decay of 56 Ni in the ejecta, 
as described by equation 4 in JET10. 
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The JET10 models show that mixing is mostly com- 
plete by the time the shock reaches the surface of the star, 
so elemental mass fractions are realistically distributed 
in radius and angle in our profiles when we spherically 
average them and map them into RAGE. We halt the 
CASTRO simulation when the shock is no less than 100 
photon mean free paths A p from the edge of the star: 

\ = — , (1) 
up 

where n is the opacity due to Thomson scattering from 
electrons, taken to be 0.288 gm cm 2 , and p is the density 
just ahead of the shock inside the star. 

This intermediate step in CASTRO allows us to cap- 
ture how heavy elements that have built up in the star 
are mixed in the initial stages of the SN without resorting 
to full 2D radiation transport models in RAGE. Such cal- 
culations, while tractable, would require too much time 
for a grid of models as large as ours and arc not nec- 
essary because photon transport does not affect mixing. 
Mixing determines the order in which emission and ab- 
sorption lines appear in the spectra over time. If metal 
lines appear soon after shock breakout in a Pop III SN 
they would be a clear signature of a core-collapse event 
because only in such explosions can mixing dredge these 
elements up from deep inside the star. Intermediate mod- 
eling in 2D is not needed for PI SN light curves because 
mixing is not vigorous enough to bring elements from 
the co re up to the photosphere of the explosion at early 
times (jJoggerst fc Whalen|[20H . 

We refer our readers to § 2 of JET 10 for further details 
on our CASTRO numerical mesh, boundary conditions, 
gravity scheme, EOS, AMR refinement criteria, and en- 
ergy deposition due to 56 Ni decay. 

2.3. RAGE 

We evolve the shock through the outer layers of the 
star, its surface, and then out into the surrounding 
medium with the rad iation hydrodynamics code RAGE 
(jGittings et al.ll2008D . RAGE (Radiation Adaptive Grid 
Eulcrian) is an AMR radiation hydrodynamics code that 
combines a second-order conservative Godunov hydro 
scheme with grey or multigroup flux-limited diffusion 
to model strongly radiating flows in ID, 2D, or 3D. 
RAGE uses atomic opacities compiled from th e exten- 
sive LANL OPLIB databasePl ljMagee et af1ll995D and can 
evolve multimaterial flows with several options for the 
EOS. We describe the physic s implemented in our RAGE 
models in iFrev et all (|2012l ) (hereafter FET12): multi- 
species advection, 2-temperature (2T) grey flux-limited 
diffusion, and energy deposition due to the radioactive 
decay of 56 Ni (| Fryer et al.ll2009fl . We note that 2T radi- 
ation transport, in which radiation and matter tempera- 
tures are evolved separately, better models shock break- 
out and its aftermath than previous IT models of SN 
explosions. We exclude gravity from our RAGE calcula- 
tions because JET10 showed that fallback onto the cen- 
tral remnant is over before breakout and that its gravity 
does not strongly affect the shock after it reaches the sur- 
face of the star. We evolve mass fractions for the same 
27 elements as in CASTRO. 

7 http:/ /aphysics2/www. t4.lanl.gov/cgi-bin/opacity/tops. pi 



2.3.1. Model Setup 

We spherically-average 2D densities, velocities, tem- 
peratures and mass fractions for the explosion, the star, 
and the circumstellar envelope from CASTRO onto a 
50,000 zone ID spherical AMR mesh in RAGE. Since we 
do not evolve radiation energy densities in CASTRO, we 
initialize them in RAGE by assuming that 

E rad = aT\ (2) 

where a is the radiation constant and T is the gas tem- 
perature in CASTRO. Also, we construct the specific in- 
ternal energy (erg gm -1 ) from T from 

3 

Egas = 2 RT i ( 3 ) 

where R = 8.314 x 10 7 erg K _1 is the ideal gas constant. 
The initial radius of the shock varies from 4.1 x 10 13 to 
1.3 x 10 14 cm in the z-series and 2.04 x 10 12 to 4.77 
x 10 12 cm in the u-series. Our root grid at setup has a 
resolution of 8.0 x 10 s cm and an outer boundary at 4.0 
x 10 14 cm for the z-series and a resolution of 1.0 x 10 9 cm 
and an outer boundary at 5.0 x 10 13 cm in the u-scrics. 
Wc allow up to five levels of refinement in both the initial 
interpolation of the profiles onto the setup grid and later 
throughout the simulation. Our setup guarantees that all 
important features in the profiles are resolved by at least 
10 grid points and can be subsampled by up to 32 times 
more points if necessary. We ensure that the explosion is 
spanned by at least 5000 zones at shock breakout so that 
the photosphere of the shock is fully resolved; failure to 
do so can lead to underestimates of luminosity during 
post-processing. 

We impose reflecting and outflow boundary conditions 
on the fluid and radiation flows at the inner and outer 
boundaries of the mesh, respectively. When the calcula- 
tion is launched, Courant times are initially short due to 
high temperatures, large velocities and small cell sizes. 
To speed up the simulation and to accommodate the ex- 
pansion of the ejecta, we periodically resample the pro- 
files onto a larger mesh as the explosion grows. Each 
rcgrid significantly increases the time step on which the 
model evolves. We remap just the explosion itself, ex- 
cluding any medium beyond the shock, and then graft 
the original environment that lies beyond this radius onto 
the shock on the new grid. We apply the same criteria 
in choosing a new root grid as in the original problem 
setup: any important density or velocity structure must 
resolved by at least 10 mesh zones and be sampled by 
up to 32 times more zones as needed, and at least 5000 
coarse grid zones are allocated to the blast after break- 
out. The outer boundary of the final largest mesh in our 
models is 5.0 x 10 16 cm, the greatest distance the ejecta 
can propagate in four months. 

2.3.2. Circumstellar Envelope 

The ambient media of z ~ 30 and z ~ 10 - 15 Pop III 
PI SNe are quite different because the former occur in 
small cosmological halos and the latter go off in primeval 
galaxies. Even low-mass Pop III st ars can photocvapo- 
rate minihalos (jWhalen et al.1 l2~008f ) . and so at z > 20 
most die in diffuse relic H II regions with n ~ 0.1 cm . 
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Fig. 1. — Kepler profiles for the shock, the star and its surrounding envelope for both u- and z-series Pop III CC SNe as initialized in 
RAGE. Left: densities. Right: velocities. 



Local densities in z ~ 15 protogalaxies are less well un- 
derstood but are likely higher. In either case, if the star 
sheds a wind it will determine the densit y profile clos- 
est to t he star. When considering PI SNe, IWhalen et al.l 
(|2012bD allowed for the possibility that massive Pop III 
progenitors have modest winds in addition to being sur- 
rounded by a diffuse H II region. However, 15 - 40 Mq 
Pop III stars are less likely to drive strong winds given 
their lower surface temperatures, luminosities, and the 
fact that stars in this mass range today manifest only 
weak winds, even at solar metallicities. Consequently, in 
this study we join a very low-mass wind profile to the 
surface of the star: 



Pw{r) 



4itr 2 vw ' 



(4) 



where to is the mass loss rate of the wind and v\y is its 
speed. The mass loss rate is just 



M t , 



tMSL ' 



(5) 



where M to t and tMSL are the total mass loss and main 
sequence lifetime of the star, respectively. We take M to t 
to be 0.01 M Q , v w to be 1000 km s" 1 , and H and He 
mass fractions to be 76% and 24%, respectively. We 
extend the wind from the surface of the star to the outer 
boundary of the mesh in RAGE. We assume that the 
wind has enough time to establish a free-streaming region 
on the grid and that any bubbles blown by the wind 
are carried far beyond the r e ach of the ejecta prior to 
the SN. As in IWhalen et al.l (|2012b[ ). for simplicity the 
wind surrounding the star is assumed to be neutral. This 
is true of stars that die as red supergiants, but not for 
blue giants, which likely ionize their envelopes. Our u- 
series light curves and spectra should therefore be taken 
as lower limits. We show initial density and velocity 
profiles for our blast models in RAGE in Figure [T] 



2.4. SPECTRUM 

To calculate a spectrum from a RAGE profile, we first 
map its densities, temperatures, velocities and mass frac- 
tions onto a 2D grid in the Los Alamos SPECTRUM 
code. SPECTRUM performs a direct sum of the lumi- 
nosity of every fluid element in this discretized profile 
to compute the total flux emitted by the ejecta along 
the line of sight at every wavelength. This procedure is 
described in FET12, and it accounts for Dopplcr shifts 
and time dilation due to the relativistic expansion of 
the ejecta. SPECTRUM also calculates the intensities 
of emission lines and the attenuation of flux along the 
line of sight, thereby capturing both limb darkening and 
absorption lines imprinted on the flux by intervening ma- 
terial in the ejecta and wind. 

Gas densities, velocities, mass fractions and radiation 
temperatures from the finest levels of refinement on the 
RAGE AMR grid are first extracted and ordered by ra- 
dius into separate files, with one variable in each file. 
Because of limitations on machine memory and time, we 
map only a subset of these points onto the SPECTRUM 
grid. We sample the RAGE radiation energy density 
profile inward from the outer boundary to determine the 
position of the radiation front, where aT 4 rises above 1.0 
erg/cm 3 . We next find the radius of the r = 20 surface 
by integrating the optical depth due to Thomson scat- 
tering inward from the outer boundary, taking kth to be 
0.288 for H and He gas at primordial composition. This is 
the greatest depth from which radiation can escape from 
the ejecta, since the minimum total opacity the photons 
would encounter is that due to just e~ scattering. 

The extracted gas densities, velocities, temperatures 
and species mass fractions are then interpolated onto a 
2D grid in r and 9 in SPECTRUM whose inner boundary 
is that of the RAGE mesh and whose outer boundary is 
10 18 cm. One hundred uniform zones in log radius are 
assigned from the center of the grid to the r = 20 surface, 
and the region between the r = 20 surface and the edge of 
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Fig. 2. — Shock breakout in the u40G (compact blue giant) run. Left panel: from left to right, velocities at 1.81 X 10 4 s, 1.87 X 10 4 s 
and 1.93 X 10 4 s. Right panel: gas temperatures at the same times. 
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Fig. 3. — Shock breakout in the z40G (red supergiant) run. Left panel: from left to right, velocities at 3.91 X 10 5 s, 3.97 X 10 5 s and 
4.14 X 10 5 s. Right panel: gas temperatures at the same times. 



the radiation front is partitioned into 5000 uniform zones 
in radius. The wind between the front and the outer edge 
of the grid is divided into one hundred uniform zones 
in log radius, for a total of 5200 radial bins. The data 
within each of these new radial bins is mass-averaged 
to guarantee that the SPECTRUM grid captures very 
sharp features from the original RAGE profile. The grid 
is uniformly divided into 160 bins in /! = cos# from -1 to 
1. Our mesh fully resolves regions of the flow from which 
photons can escape the ejecta and only lightly samples 
those from which they cannot. 



The sum of the luminosities over all wavelengths in one 
spectrum yields the bolometric luminosity of the SN at 
that moment. Many such luminosities computed over a 
range of times constitutes the light curve of the explo- 
sion. We cover shock breakout with 50 spectra uniformly 
spaced in time and the rest of the light curve with 200 
spectra that are logarithmically distributed in time out 
to 4 months. Each SPECTRUM calculation requires 3 - 
6 hours on 32 processors on LANL platforms. 

3. BLAST PROFILES, LIGHT CURVES AND SPECTRA 
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Fig. 4. — The freely expanding u40G SN at intermediate times, 
and 8.74 X 10 5 s. Right panel: densities at the same times. 

Hydrodynamical profiles for shock breakout are shown 
for the fiducial cases u40G, a blue giant progenitor, and 
z40G, a red supergiant progenitor, in Figures [2] and [3] 
The reverse shock in the z40G velocity profile in Figure 
[U which vigorously mixe s the interior of the star prior 
to breakout (Figure 5 of iJoggerst et al.l I2O1O0 . is gone 
by the time the shock reaches the surface of the star at 
3.91 x 10 5 s. With no reverse shock there is no fur- 
ther Rayleigh- Taylor mixing, and the metals are essen- 
tially frozen in mass coordinate as the ejecta expands. 
Our procedure for spherically averaging 2D CASTRO 
mass fractions to approximate their radial distribution 
in ID in RAGE therefore captures all mixing in the ex- 
plosion. The SN is not visible to an external observer 
before breakout because electrons between the surface 
and the shock scatter the photons. 

When the shock breaks through the surface of the star 
its velocity abruptly doubles (z40G) or triples (u40G). 
As photons stream free from the shock they blow off the 
outer layers of the star. The breakout transient is visible 
as the flat plateau in gas temperature that moves ahead 
of the ejecta. This plateau traces the leading edge of the 
radiation front. The height of the plateau is the tem- 
perature to which the radiation pulse heats the gas as 
it passes through it, not the temperature of the fireball 
itself, which is much higher. The u40G shock is much 
hotter than the z40G shock at breakout (~ 100 eV vs. 
10 eV) because it has broken out of a much more com- 
pact star and done less PdV work on its surroundings. 
Consequently, the radiation front initially heats the wind 
to 20 eV in u40G but to only 1 eV in z40G. As the shock 
expands its spectrum softens, and the temperature of the 
gas behind the radiation front falls over time: from 20 
eV to 8 eV in u40G from 1.81 x 10 4 s to 1.93 x 10 4 s 
and from 1.2 eV to 0.4 eV in z40G from 3.91 x 10 5 s to 
4.14 x 10 5 s. 

In both explosions, shock breakout coincides with ra- 
diation breakout because the moderate wind density at 
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panel: from left to right, velocities at 4.96 X 10 5 s, 6.80 X 10 s s 

the surface of the star cannot confine the radiation front. 
As both shocks descend the r -2 wind gradient, the u40G 
SN accelerates but the z40G SN speeds up and then de- 
celerates. The u40G shock continues to accelerate even 
though ambient densities are 1000 times higher at the 
surface of the star because the explosion is much hotter 
and the shock is more radiatively driven than the much 
cooler z40G shock. Soon after breakout, both SNe as- 
sume self-similar free expansion velocity profiles and ho- 
mologously expand, as we show in Figure U for u40G. 
This self-similar expansion continues until the end of 
the simulations at 4 months. Eventually, spherical di- 
lution of the ejecta renders the SN transparent, and it 
enters the nebular phase. At this stage, fluid elements 
in the ejecta may no longer remain in radiation equilib- 
rium with their neighbors so the true opacity may deviate 
from the OPLIB opacities, which assume local thermo- 
dynamic equilibrium (P TE). Also, the K irchhoff-Planck 
relation (equation 6 of iFrev et al.l I2012T ) may not fully 
hold in the SPECTRUM code. There may therefore be 
some inaccuracy in our spectra at late times, but these 
events are only visible in the NIR at much earlier stages 
of the explosion, as we later show. 

3.1. Light Curves / Spectra 

We show bolometric luminosities for our Pop III CC 
SNe in Figure [5j Peak luminosities at shock breakout 
vary from 8 x 10 44 to 1.5 x 10 46 erg s" 1 in the z-series 
and from 3 x 10 44 to 5 x 10 45 erg s _1 in the u-series. 
At early times the light curve is powered primarily by 
the conversion of kinetic energy into thermal energy by 
the shock, so for a given progenitor mass the peak lu- 
minosities increase with explosion energy. They peak at 
earlier times with greater E because the shock reaches 
the surface of the star in less time. Breakout also gener- 
ally happens later with more massive stars because they 
have greater radii in both series, with the exception of 
the 25 Mq star because it dies with the smallest radius, 
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Fig. 5. — Bolometric light curves for 17 Pop III CC SNe in the source frame. Upper left: all eight z-series light curves (red supergiant 
progenitors). Upper right, lower left and lower right are u!5, u25, and u40 (blue giant progenitor) light curves, respectively. 



as shown in Figure [T] For a given progenitor mass and 
E, u-series SNe are somewhat less luminous than z-series 
SNe, and this is due to a tradeoff between shock temper- 
ature at breakout and the radius of the star: 

L = 47rer 2 crT 4 . (6) 

Here, a is the Stefan-Boltzmann constant, e is the grey- 
body correction to the blackbody luminosity assumed for 
the shock (e ~ 0.1), and T is the temperature of the shock 
at the TTh = 1 surface, where kth = 0.288. The red z- 
series stars have radii 10 - 30 times greater than u-series 
stars of equal mass but their shocks have lower tempera- 
tures at breakout because they must do more PdV work 
before crashing through the surface of the star. Their 
respective temperatures are evident in their spectra at 
breakout, which we show in Figure [5] The u40G model 
has much more spectral energy below 10 A than z40G at 



breakout because the fireball has a much higher tempera- 
ture. The duration of the breakout transient is governed 
by the light-crossing time of the star and is 10 - 30 times 
greater in the z-series than in the u-series because of their 
larger radii. 

As the shock expands it cools, and its emission at later 
times is powered by the radioactive decay of 56 Ni. The 
duration of this emission can be approximated by the 
radiation diffusion timescales in the ejecta: 

t d ^ K^M^E-i, (7) 

where M e j is the mass of the ejecta, K is the average 
opacity of the ejecta, and E is the explosion energy. The 
u-series SNe dim after about 3 months but the z-series 
explosions exhibit much more extended emission reminis- 
cent of Type Hp SNe, whose progenitors are also thought 
to be red supcrgiants. In particular, the z40 series remain 
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Fig. 6.— Spectral evolution of the fireball. Right: z40G at 3.92 X 10 5 s (black), 5.55 X 10 s s (blue), 1.49 X 10 6 s (red), 8.66 X 10 6 s 
(green) and 2.95 X 10 7 s (purple). Left: u40G at 1.87 X 10 4 s (black), 6.82 X 10 4 s (blue), 2.06 X 10 s s (red), 1.21 X 10 6 s (green) and 
9.69 X 10 6 s (purple). 



above 10 42 erg s _1 for up to 300 days. However, because 
the shock cools over time, the region of its spectrum that 
is redshifted into the NIR in the observer frame dims 
below detection limits well before 1 year, as we discuss 
below. 

We show the evolution of the z40G and u40G spectra 
in Figure [51 Unlike with much more energetic Pop III 
PI SNe, the outer regions of the ejecta and the envelope 
are never fully ionized, so they always blanket the spec- 
trum with absorption lines at the short wavelength limit. 
As the shock expands and cools its surrounding envelope 
imposes more lines on the spectrum and continuum emis- 
sion at long wavelengths slowly rises. The flux at longer 
wavelengths increases with time because the ejecta cools 
and its surface area grows. From these profiles, it is 
clear that fitting a blackbody profile to a bolometric lu- 
minosity to approximate a spectrum can lead to serious 
overestimates of flux at short wavelengths, because it er- 
roneously reinstates luminosity that is actually removed 
by line blanketing. This caveat is especially pertinent to 
detection thresholds for high-redshift SNe because flux 
from this region of the spectrum is a principal component 
of the NIR signal of the explosion in the observer frame. 
The collision of the shock with a realistic circumstellar 
envelope (which sets its temperature) and the opacity of 
the ejecta and envelope also crucially shape the spectra 
in the source frame, and thus the NIR light curve in the 
observer frame. 

4. POP III CC SN DETECTION THRESHOLDS 

We calculate NIR light curve s for our SNe w ith the 
photometry code developed by ISu et al.l (|2011l) . Each 
spectrum is redshifted before removing the flux that is 
absorbed by intervening neutral hy drogen a l ong th e line 
of sight according to the method of iMadaul (|1995f ). The 
spectrum is then dimmed by the required cosmological 
factors. We linearly interpolate the least sampled data 
between the input spectrum and filter curve to model the 
light curve in a given filter. 



4.1. NIR Light Curves 

At every redshift for each supernova, we calculate the 
NIR signal in JWST NIRCam filters above and below 
the Lyman limi t to find the filter in which the SN is 
brightest. As in iWhalen et al.l (|2012d ). we find that the 
CC SNe in our study are most luminous just redward of 
1216 A in the source frame at the redshifts we consider. 
We show NIR light curves for z-series and u-series SNe 
at z = 7, 10 and 15 in Figures and U The NIRCam 
photometry limit is AB magnitude 32, and so six of the 
eight z-series SNe will be visible for 20 - 600 days at z = 
7 but only one of the u-series explosions, u40G, will be 
visible, and only for ~ 30 days. Five of the z-series SNe 
will be visible at z = 10: zl5B for less than 1 day, zl5D 
for 75 days, zl5G for 220 days, z40B for 275 days and 
z40G for over 400 days. At z = 15 the zl5G SN is visible 
for 250 days and the z40B and z40G explosions can be 
seen for 300 days. The peak in emission for all the SNe 
advances to later times and is broader at earlier epochs 
because of cosmological rcdshifting. For a given mass, 
the NIR luminosity increases with the explosion energy. 

Shock breakout is not visible from earth in any of these 
SNe because the x-rays and hard UV in the transient are 
absorbed by neutral H at high redshift. The flux that 
is redshifted into the NIR varies much more rapidly (on 
timescales of ~ 100 days) than the bolometric luminosity 
(3-7 years) in the observer frame because the spectra 
evolve as the fireball expands. The light curves rise more 
rapidly than they fall, so they are most recognizable as 
transients in their early stages. However, given JWST 
survey times of 1 - 5 years, it will be possible to identify 
these events as SNe at any stage above detection thresh- 
old. Because they are much dimmer than Pop III PI SNe 
and because z ~ 7 galaxies will be more luminous than 
z ~ 10 - 15 protogalaxies, some Pop III CC SNe may be 
somewhat more challenging to discriminate from their 
host galaxies (with which they likely overlap in color- 
color space). However, they will still be more luminous 
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Fig. 7. — NIR light curves for 2 = 7 Pop III CC SNe at 1.63 fim in the observer frame. Left: z-series. Right: u-series 



than most galaxies of that epoch and may be much more 
frequent than PI SNe if baryons in some metal-free halos 
collapse and fragment into small swarms of stars rather 
than a few very massive ones, as the latest numerical 
simulations suggest. 

The prospects for detecting Pop III CC SNe with 
JWST at z > 10 are better than those for PI SNe (e.g., 
iHummel et al.l 120121) because of their higher rates, but 
they will be too dim to be found in future all-sky NIR 
surveys such as Euclid, whose detection limit is AB mag- 
nitude 24 at 2 fim, or WFIRST and WISH, whose pho- 
tometry limits will be AB magnitude 27 at 2.2 /xm. They 
also will fall below the detection threshold in the Y-band 
centered at 1.0 /mi for the Large Synoptic Survey Tele- 
scope (LSST) and the Panoramic Survey Telescope & 
Rapid Response System (Pan-STARRS), whose photom- 
etry limits will be at most AB magnitude 25 and 27, 
respectively. Our SPECTRUM calculations show that 
these SNe will be even dimmer in the optical in LSST and 
Pan-STARRS because of extinction by the IGM at high 
rcdshift. We note that while we have considered only Pop 
III CC SNe, our detection limits hold for progenitors of 
any metallicity. The central engine of the explosion de- 
pends mostly on the entropy profile of the inner 3-4 Mq 
of t he star, which does not vary strongly with met allic- 
ity (IChieffi fc Limongil[200l IWooslev fe Hegerll2007l) (see 
also Figure 1 in lWhalen fc Frverj^OiaT The light curves 
and spectra of CC SNe at higher metallicities should 
therefore be bracketed by those of our compact blue giant 
and red supcrgiant Pop III progenitor stars. 

5. CONCLUSION 

Although direct detections of Pop III CC SNe will not 
probe the very first stellar populations, they will reveal 
the properties of stars in z ^ 10 - 15 galaxies and trace 
their formation ratesand early galactic chemical evolu- 
tion in general. Furthermore, because they will likely be 
brighter than most primitive galaxies of that era, explo- 
sions of low-mass Pop III stars (and indeed any CC ex- 
plosion) will reveal the existence of primeval galaxies on 
the sky that might otherwise escape detection by JWST 
or future ground-based 30-metcr class telescopes. We 



note that while we have considered explosion energies of 
0.6 - 2.4 Be, higher energies may be possible for some 
CC SNe and would be visible at even higher redshifts. 

There are also scenarios in which CC explosions can 
produce luminosities that rival those of PI explosions 
that could be visible at much earlier epochs. For exam- 
ple, if the shock collides with a dense shell ejected during 
a luminous blue variable (LBV) outburst prior to the SN, 
an extremely bright event in the UV can result that could 
be detected above z ~ 15 (e.g. ISmith fc McCravl 120071; 
Smith et all 120071; Fan Marie et all 120101: IMoriva et al.l 



20101 , 12012t iTanaka et all 12012ft . We are now studying 



the observational signatures of these su perluminous Pop 
III Type Iln SNe (jWhalen et al.ll2012aft . There may also 
be hypernovae, very energetic explosions of 40 - 60 Mq 
stars with energies of 10 erg that are intermediate to 
those of CC and PI SNe. This kind of explosion, whose 
existence is inferred in part from the elemental abun- 
dances imprint ed on a few of the m ost metal-poor stars 
found to date (|Iwamoto et al.|[2~005f ). could be detected 
at redshifts that bridge those at which CC and PI SNe 
can be found, 15 < z < 20. We are also calculating light 
curves and spectra for such events. 

Strong gravitational lensing by massive galaxies and 
clusters at z ~ - 1 could magnify flux from Pop 
III supernovae, impr oving prospects for their detection 
(|Rvdberg et al.ll2012t) . The probability that flux from a 
Pop III SN would be boosted in an all-sky survey and 
its magnification depend on the event rate on the entire 
sky at a given rcdshift. We have performed preliminary 
Markov Chain Monte Carlo ray-tracing calculations that 
suggest that the probability that a z ~ 20 event is lensed 
is ~ 1 - 5% for flux boosts of 2 - 5. Much higher magnifi- 
cations (10 - 300) arc possible near the edges of massive 
clusters but the search volumes and probabilities of en- 
countering high-z SNe are much smaller. We continue to 
study strong lensing of z ~ 20 events, the highest red- 
shifts ever attempted, in order to assess its potential to 
discover primeval SNe and galaxies. 

Can later stages of the SN remnant be detected? 
IWhalen et"al~l (| 20081) found that most of the energy of 
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15 and 40 M© Pop III SNe is eventually radiated away 
as H and He lines as the remnant sweeps up and shocks 
pristine gas. At lower redshifts this energy would in- 
stead be lost to fine structure cooling by metals. In cither 
case, the emission is too diffuse, rcdshiftcd and drawn out 
over time to be detected by any upcoming instruments. 
Also, unlike PI SNe, CC SNe deposit little of their energy 
into the cosmic microwave backgro und (CMB) by inverse 
Compton scattering at z ~ 20 (jKitavama fc Yoshidal 
120051: iWhalen et all 120081) . and even less at lower red- 
shifts because the density of CMB photons falls with cos- 
mological expansion. Consequently, early populations of 
low-mass Pop III SNe will probabl y not imp ose excess 
power on the CMB at small scales (|Oh et al.ll2003l ). For 
the same reason it will probably not be possible to di- 
rectly image Sunyaev-Zeldovich fluctuations from indi- 
vidual Pop III CC SN remnants with the Atacama Cos- 
mology Telescope or the South Pole Telescope. 

Although their event rates make it unlikely that Pop 
III CC SNe will be detected in absorption at 21 cm 
at z > 10, new calculations reveal that enough syn- 
chrotron emission from their remnants is rcdshiftcd into 
the 21 cm band above z ~ 10 to be directly detected by 
the S quare Kilometer Array (SKA) ([Meiksin &: Whaler] 
120121 ). Somewhat more energetic hypernovae could be 
detected by existing facilities such as the Extended Very- 
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